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Abstract 

We investigate degenerate saddle point problems, which can be viewed as limit cases of standard mixed formulations of symmetric 
problems with large jumps in coefficients. We prove that they are well-posed in a standard norm despite the degeneracy. By 
wellposedness we mean a stable dependence of the solution on the right-hand side. A known approach of splitting the saddle 
point problem into separate equations for the primary unknown and for the Lagrange multiplier is used. We revisit the traditional 
Ladygenskaya-Babuska-Brezzi (LBB) or inf-sup condition as well as the standard coercivity condition, and analyze how they 
are affected by the degeneracy of the corresponding bilinear forms. We suggest and discuss generalized conditions that cover the 
degenerate case. The LBB or inf-sup condition is necessary and sufficient for wellposedness of the problem with respect to the 
Lagrange multiplier under some assumptions. The generalized coercivity condition is necessary and sufficient for wellposedness of 
the problem with respect to the primary unknown under some other assumptions. We connect the generalized coercivity condition 
to the positiveness of the minimum gap of relevant subspaces, and propose several equivalent expressions for the minimum gap. Our 
results provide a foundation for research on uniform wellposedness of mixed formulations of symmetric problems with large jumps 
in coefficients in a standard norm, independent of the jumps. Such problems appear, e.g., in numerical simulations of composite 
materials made of components with contrasting properties. 

Key words: Wellposedness, mixed, symmetric, saddle point, Lagrange multiplier, Ladygenskaya— Babuska-Brezzi (LBB) condition, inf-sup 
condition, coercivity, minimum gap between subspaces. 
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1. Introduction 

Degenerate saddle point problems, e.g., can be viewed as 
limit cases of mixed formulations of symmetric problems 
with large jumps in coefficients, corresponding to an infi- 
nite jump. We prove that the degeneracy docs not affect the 
wellposedness in a standard norm under some natural as- 
sumptions, using ideas that arc initiated by [3, 4, 5, 6, 7, 14, 
15]. By wellposedness, contrary to illposcdness, we mean a 
stable dependence of the solution on the right-hand side. 
Results of this paper provide a foundation for research on 
uniform wellposedness of mixed formulations of symmet- 
ric problems with large jumps in coefficients in a standard 
norm, independent of the jumps. 
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The necessary and sufficient condition, e.g., [9, 10], of 
the standard wellposedness of an operator equation with 
an arbitrary right-hand side is the existence of a bounded 
inverse of the operator. We argue that in some practical 
cases the equation is degenerate, i.e. the inverse operator 
does not exist. Assuming that the right-hand side is in the 
operator range, a solution exists, but is not unique. To make 
the solution unique we factor out the operator null-space. 
This leads to a natural generalization, where boundedness 
of the pseudoinverse of the operator is used as the necessary 
and sufficient condition of wellposedness of a degenerate 
operator equation, by analogy with [13, 15]. 

With this idea in mind, we revisit necessary and sufficient 
conditions of wellposedness of an abstract mixed problem. 
In the symmetric case we consider here, the mixed problem 
can be interpreted as a variational saddle point problem. 
For generalized saddle point problems we refer the reader, 
e.g., to [11]. 

We start in Section 2 with a standard abstract symmetric 
mixed problem as in [9, 10]. By analogy with [14, 17], we 



split the saddle point problem into two equations, for the 
primary unknown and for the Lagrange multiplier. This 
split is somewhat implicit in [9, 10]. The equation for the 
primary unknown is self-consistent, since here we eliminate 
the Lagrange multiplier from the mixed system using an 
orthogonal projector. 

Following, e.g., [10], we discuss the traditional neces- 
sary and sufficient conditions of wellposedness, namely, the 
Ladygenskaya-Babuska-Brezzi (LBB) or inf-sup condition 
and the coercivity condition. The LBB or inf-sup condition, 
considered in Section 3, is necessary and sufficient for a sta- 
ble dependence of the Lagrange multiplier on an arbitrary 
right-hand side. 

We review the traditional point of view that the coerciv- 
ity condition is a necessary and sufficient condition of well- 
posedness of the problem. In Section 4, an operator form 
of the dual variational problem without assuming the co- 
ercivity condition is considered. We examine the unique- 
ness of the solution and describe all possible multiple solu- 
tions for a given right-hand side. All admissible right-hand 
sides are determined. We formulate several equivalent nec- 
essary and sufficient conditions of wellposedness in terms 
of closedness of relevant subspaces. We also derive a geo- 
metrical condition — a positiveness of a minimum gap [12] 
between relevant closed subspaces. 

A possible application of our theory is the Hcllingcr- 
Rcissncr formulation, e.g., [1], of nonhomogeneous Lame 
equations for media with (almost) rigid inclusions, where 
the Lagrange multiplier is the displacement, and we get an 
operator equation for the stress on the closed subspace of 
divergence free (in a weak sense) stresses. Infinitely large 
Lame coefficients A and /x, in a subdomain, result in a 
null-space of the operator in the equation for the stress, 
so the inverse operator does not exist and the problem is 
not wcllposcd in a traditional sense. Our abstract geomet- 
rical condition of generalized wellposedness in this example 
is equivalent to a possibility of extension of displacements 
preserving the energy norm of the Lame operator. It has 
been proved in [4, 7] that such an extension is possible un- 
der some assumptions. We expect that in the limit case of 
infinitely large Lame coefficients A and [i in a subdomain 
the pseudoinverse of the operator is bounded, which makes 
the problem wellposed for the stresses in the sense, i.e. 
the L2 norm of the stress is stable even if the Lame coef- 
ficients are large in a subdomain. We plan to address this 
application in the future. 

2. Abstract symmetric saddle point problems 

In this section we essentially follow well known argu- 
ments, e.g., [10], with some simplifications due to the sym- 
metry of the saddle point problem and our unwillingness 
to introduce dual spaces. Straightforward manipulations, 
using a pair of complementary closed subspaces, allow us, 
as in [14, 17], to formulate separate equations for the pri- 
mary unknown and for the Lagrange multiplier of the sad- 



dle point problem; sec, e.g., survey [8, Sec. 6] for similar 
matrix null-space methods. We start by formulating and 
investigating the problem using bilinear forms, and then 
repeat the arguments for operator-based formulations that 
are used in the last section of the paper. 



2.1. Formulations using bilinear forms 

Let H and V be two real Hilbert spaces with scalar prod- 
ucts and norms denoted by (•, -)h, || ■ ||h and (•, -)v, || ■ ||v 
correspondingly. Let a(-, •) : H x H — > R and &(■, ■) : 
H x V — > R be two continuous bilinear forms with a(-, ■) 
symmetric and nonncgative definite. We consider the fol- 
lowing problem: for a given g G H and / G H find a G H, 
called the "primary unknown," and u G V, called the "La- 
grange multiplier," such that 

a(a,e) + b{e,u) = (g,e) n , Ve G H, 
b(a -f,v) = 0, V?; G V. 

We place the right-hand side / "inside" of the form b as it 
allows us to take / G H, not to introduce the dual space 
V, and makes several statements somewhat simpler. We 
call (1) a saddle point problem, since equations (1) are the 
optimality conditions and their solution is a saddle point 
for the Lagrangian, e.g., [10], defined by a(cr,a) + 2b(a — 

f, u ) ~ 2(g,cr)H- 

We call a linear manifold, not necessarily closed, a "sub- 
space" and a closed linear manifold a "closed subspace." 
Let us introduce a special notation N C H for the closed 
subspace, which is the null-space of the bilinear form &(•, •) 
with respect to its first argument, i.e. N = {e G H : 
b(e,v) = 0, Vw G V}. Let us denote by P = N 1 C H the 
closed subspace which is H-orthogonal (complementary) 
to N. Closed subspaces N and P play important roles in 
this paper, so let us introduce an H- orthogonal projector P 
on H such that N(P) = N and R(P) = P and the com- 
plementary projector P = I — P with R(P- L ) = N and 
N(P ± ) = P, where by R(P) we denote the range of opera- 
tor P and, with a slight abuse of the notation, by N(P) we 
denote the null-space of operator P. We assume through- 
out the paper, unless stated otherwise, that a bounded op- 
erator is defined everywhere on a corresponding space. As 
an orthogonal projector, operator P : H — > H is bounded 
H-selfadjoint, P = P*, and satisfies P = P 2 . 

In the first equation of system (1), let us split it into two 
equations, by plugging separately e = Pe G P and e = 
P^e G N and using the fact that 6(P _L e, u) = 0, Ve G H. 
The second equation in system (1) has a simple equivalent 
geometric interpretation: a — f G N, or [a — f, e)n = 
0, Ve G P. We then rewrite system (1) in the following 
equivalent form: 



a(a, e) + b(e, u) = (g, e) H , Ve G P, 

a(a,e) = (g,e)n, Ve G N, (2) 
(cr-/,e) H =0, VeeP. 

Now we make an important observation that we can treat 
the first line in system (2) as an equation for the Lagrange 
multiplier u, given the primary unknown cr, i.e. 

K e > u ) = (.9, e)h — a(a,e), Ve G P. (3) 

The last two lines in system (2) involve neither the Lagrange 
multiplier u, nor the bilinear form b, and can be used to 
determine the primary unknown cr: 

a(a,e) = (g,e) H , Ve G N, 
(cr-/,e) H = 0, VeeP. 

System (4) describes, e.g., [10], the optimality conditions 
of the constrained minimization problem inf {a(cr, cr) — 
2(3, <7)h}, a G H : (d - /, e) H = 0, Ve G P. 

2.2. Operator-based formulations 

In addition to the formulations above involving bilin- 
ear forms, it is convenient to consider equivalent operator- 
based formulations. We associate with the forms a and b two 
linear continuous operators A : H — ► H and B : H — > V 
defined by (Aer, e)n = a(cr, e), (Ba,v)v = b(a,v), Ve, cr G 
H, f G V. In this definition of A and B we follow a 
slightly simplified, e.g., [11, 17], rather than standard [10], 
approach, namely, we do not need dual spaces H' and V . 
Now, we reformulate the main statements of subsection 2.1 
using the just defined operators A and B. The following 
operator formulation 

!Aa + B*u = q in H, 
(5) 
B{a - /) = in V 

is equivalent to the original problem (1) with the bilinear 
forms, where the adjoint operator B* : V — > H is defined, 
as usual, by (cr, _B*u)h = {Ba,v)v, Vcr G H, v G V. The 
operator A is selfadjoint and nonnegative definite, A = 
A* > on H since it is defined by the symmetric and 
nonnegative definite form a. 

We notice that the second equation in system (5) has 
the same geometric interpretation as in the case of bilinear 
forms-based system (1): cr — / G N(£?). The null-space 
N(B) C H and its H-orthogonal complement R(B*) C H 
have already been denoted by N and P, correspondingly, 
and introduced together with the H- orthogonal projector 
P on H such that N = N(P) = N(B) and P = R(P) = 
R(P*) in the previous subsection. 

We split the first equation in system (5) in two orthogonal 
parts corresponding to N and P, using that PB*u = B*u 
and P^B*u = 0, since R(B*) C P. We replace B with P, 
since they share the same null-space, in the second equation 



in system (5) to get the following equivalent form of system 
(2): 

{PAo + B*u = Pg in H, 
P x Aa = P^g in H, (6) 
P(cr - /) = in H. 

We notice that the first line in system (6) is an equation 
for the Lagrange multiplier u, given the primary unknown 
cr, as in (3), i.e. B*u = P(g — Aa). 

We next discuss the necessary and sufficient conditions 
from [10] of wcllposedness of the problem and make it clear 
why one can find weaker necessary and sufficient conditions. 
To simplify our arguments, we take advantage in the rest 
of the paper of the split of the original system into separate 
equations for the Lagrange multiplier u and the primary 
unknown cr that we have described in this section. It is 
important to realize, however, that we have not made any 
substitutions, neither in the solutions u and cr, nor in the 
right-hand sides / and g. So whatever statements we next 
prove concerning the dependence of the solutions u and 
cr on the right-hand sides / and g, these statements are 
equally applicable to both the separate equations and to the 
original system in either bilinear form- or operator-based 
context. 



3. Inf-sup or LBB condition 

In this section, we discuss a traditional assumption, be- 
ing recently referred to as Ladygenskaya-Babuska-Brezzi 
(LBB) condition, see Babuska and Aziz [2], Brezzi and 
Fortin [10], Ladyzhenskaya [16], that the range of operator 
B : H — > V, denoted by R(B), is closed. The closedness of 
a range of a closed operator is ultimately connected to the 
boundedness of the operator (pseudo-)inverse, e.g., [12]. 

In our specific situation, operator B is bounded with 
the closed domain H and, thus, is closed, so its (pseudo- 
inverse B^ 1 : R(-B) — > H/N(.B) is also closed. It is neces- 
sary to use a factor-space here to define the inverse, since 
the standard operator inverse B~ x : R(£>) — > H does not 
exist if N(_B) is nontrivial. We note that N(£?) is closed 
and that the factor-space H/N(_B) is a Hilbcrt space, as 
is H. In a Hilbcrt space, a convenient set of representants 
for the classes in the factor-space is simply the correspond- 
ing orthogonal complement, e.g., H/N(£?) is isomctrically 
isomorphic to P = (N(B)) ± C H, so we set ||cr|| H / N ( B ) = 
II Per || h- The subspace R(B) is the domain of the closed op- 
erator B- 1 : R(S) -> H/N(S) therefore, R(B) is closed 
if and only if B^ 1 : R(B) H/N(P) is bounded. Closed- 
ness of R(B) is equivalent to closedness of R(B*), so all 
the arguments above can be equivalcntly reformulated for 
the adjoint operator B* and its (pseudo-)inverse. 

When written in terms of inequalities involving the bi- 
linear form b : 



inf sup 



b(o~, v) 



inf 



\Ba\ 
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» eH »ev IM|h/n(b)IMIv ct 6 h ||ct||h/n(b) 

1 



or, equivalently, 



b(a,v) 



inf sup - 

»eV ff£H ||CT||H||W||V/N(i 



\B~ l \ 



= inf 



> 0, 



R(B)-»H/N(B) 



H 



»6V F||V/N(B* 
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R(B*)^V/N(B*) 



>o, 



the LBB condition is also known as the inf-sup condition, 
see Babuska and Aziz [2], Brezzi and Fortin [10], where 
V/N(S*) means the factor-space of V with respect to the 
closed subspace N(B*). We implicitly assume that the ar- 
guments in the inf-sup formulas above and throughout the 
paper do not make both the numerator and the denomi- 
nator vanish. In Ladyzhenskaya [16], the inf-sup condition 
does not appear to be explicitly formulated, instead, closed- 
ness of a range of the gradient operator is investigated in 
connection with wellposcdncss of the diffusion equation. 

We note that the induced norms of an operator and its 
adjoint are equal, so both inf-sup expressions above are 
equal to the same constant that we call q,. If at least one 
of the spaces H or V is finite dimensional then the value Cf, 
is positive automatically, so it becomes important how Cb 
depends on some parameters, e.g., on the dimension. 

Let us mention that in many practical applications the 
space V can be naturally defined such that N(£?*) = {0}, 
so the latter inf-sup expression of the LBB condition takes 
the form 



b{a,v) 



inf sup 
"ev^gH ||ct||h|M|v 



c b > 0, 



which can be most often seen in publications on the subject. 
We now contribute our own equivalent formulations of the 
LBB condition. 

Lemma 3.1 Subspaces R(B) C V and R(BB*) C V are 

closed simultaneously. Moreover, if either of them is closed 
we have R(BB*) = R(B). 

Proof. If BB*v = then (B*v, B*v)n = 0, i.e. B*v = 0, 
which proves that N(BB*) = N( B*). Taki ng an ort hogo- 
nal complement to both parts gives R(BB* ) = R(-B) as the 
operator BB* is selfadjoint. Triviall y, R(B B *) C R(B ). 
If the range R(BB*) is closed then R(B) = R(BB*) = 
R(BB*) C R(B), but clearly R(B) C R(B), which proves 
closedness of R(B) = R(BB*). 

To prove the inverse statement, assuming that R(B) 
is closed, we invoke the orthogonal decomposition argu- 
ment 2 : H = R(B*) © (R(B*)) 1 - = R(B*) © N(B) since 
R{B) and thus R(B*) are closed. Multiplying this equal- 
ity by B gives R{B) = BH = B(R(B*) © N(5)) = 
BR(B*) = R(BB*). □ 



We use the previous lemma to introduce (BB*) 1 : 
R(BB*) -> V/N(B*) in the next Lemma 3.2. It is nec- 
essary to use the factor-space V/N(£?*) here, since the 
standard inverse (BB*)^ 1 : R(BB*) — > V does not exist 
if N(B*) is nontrivial. 

Lemma 3.2 Closedness of R{B) C V is equivalent 
to boundedness of the operator (BB*)^ 1 : R(BB*) — ► 
V/N(B*). 

Proof. By Lemma 3.1, closedness of R(B) C V is equiv- 
alent to closedness of R(BB*) C V. We use several well- 
known statements on closed operators, e.g., [12], applied 
to the operator BB*, that we have already reviewed in the 
second paragraph of this section for the operator B. The 
operator BB* is bounded and has the closed domain V, so 
the operator is closed and its (pseudo-) inverse (BB*)^ 1 : 
R(BB*) -> V/N(B*) with the domain R(BB*) C Vis 
also closed. The domain R(BB*) C V of the closed opera- 
tor B- 1 : R(BB*) -> H/N(B) is closed if and only if the 
operator is bounded. □ 

If R(B) is closed then, using Lemmas 3.1 and 3.2, 
R(B) = R(BB*) and we can derive the following useful 
formula 



P = B*{BB*y 1 B : H 



H. 



(7) 



This proof is suggested by an anonymous referee 



Indeed, we first note that R^BB*)- 1 ) C V/N(B*) is 
multiplied by B* in (7), so the product is independent 
of the choice of a representant from the equivalence class 
V/N(B*) and, thus, is correctly defined. Second, righ-hand 
side of (7) is a linear and bounded operator as a product of 
linear and bounded operators. Moreover, it is an orthogonal 
projector on H since it is selfadjoint and idempotent, and 
has the null-space the same as the orthoprojector P has. 

If the LBB condition is not satisfied, i.e. R(B) is not 
closed, then the domain of definition of the operator 
B*(BB*y 1 B is the subspace R(B*) © N(B), which is 
not closed, and formula (7), where P is the orthogonal 
projector on H with N(P) = N(B), clearly does not hold. 

Let us note that in the case of finite dimensional spaces 
H and V the range R(B) is evidently closed, the opera- 
tor (£?*)+ = (BB*) + B is the well-known Moore-Penrose 
pseudo inverse of B* , and P = B*(B*) + is the well known 
formula for the orthogonal projector onto the range of B* . 

If a is an exact solution of system (5), then u in (5) can be 
found from the equation B*u = —Aa+g E R(B*). If a is an 
approximate solution of system (5) such that the condition 
Aa — g G R(B*), which is necessary and sufficient for the 
existence of u, does not hold, then u can be computed from 
the projected equation B*u = P(—Aa + g) G P. Both the 
original and the projected equations for u are wellposcd by 
the LBB assumption, i.e. R(B*) = P and 

Hill , 1 ii 

U v/N(B-) < — k H + - h H- 

c b c b 

Whether the LBB assumption is necessary for wcllposed- 
ness of the equation for u depends on if the set of all possible 
right-hand sides g — Aa gives the whole subspace R(B*), 



i 



see [10]. For example, in a practically important case g = 
we have B*u = -Ao = -PAr G R{PA) C R(P). If the 
latter inclusion is strict, it opens up an opportunity for a 
weaker, compared to the original LBB, assumption of well- 
poscdness of the above equation for u. 

In the present paper, however, we are concerned with 
finding a, not u. The LBB condition for the bilinear form 
b appears to be of no importance for our results in the 
next section where we analyze wellposedness of system (5) 
with respect to the a unknown only, assuming that the u 
unknown is of no interest, or can be found for a given a 
using some postprocessing. 

4. Coercivity conditions 

4.1. The standard coercivity condition 

We finally get to the main topic of the paper: an assump- 
tion on A which is a condition of wellposedness of (5) with 
respect to a. For the reader's convenience, we briefly repeat 
the necessary notation and the system of equations for a 
to make this section self-consistent. Let H be a real Hubert 
space and P be an orthoprojector in H with a null-space 
N(P) = N and a range R(P) = P — we emphasize that the 
range of any orthoprojector in a Hilbert space is closed. Let 
A be a linear and bounded operator such that < A* = A 
on H. The last two lines in system (6) represent an oper- 
ator form of system (4) ; they do not involve the Lagrange 
multiplier u or the operator B and determine the primary 
unknown a G H: 

{ P ± (Act - g) = in H, 

(8) 

{ P(<r - /) = in H, 

where jgH and / G H are given and = I — P. We can 
also replace system (8) with the following equivalent single 
equation: 

P^A | N i> = P^g - P^APf G N, a = tp + Pf, (9) 

where in (9) we take a restriction of the operator P L A on 
its invariant closed subspacc N, and we are looking for a 
solution ip G N. Then the necessary and sufficient condition 
of wellposedness of problem (9) for an arbitrary g G H is, 
clearly, that the range of P^A |n is N. This leads to the 
traditional assumption, see [10], a(a,a) > c a > 0,Vcr G 
N, ||c||h = 1 or, in an operator form, A > c a I on N C H, 
since A is selfadjoint nonnegative. Thus, this assumption 
is also necessary and sufficient [9, 10] for wellposedness of 
system (5) with respect to a for an arbitrary g G H. In 
the rest of the section, we analyze the scenario, where A 
is selfadjoint nonnegative on H, but may be degenerate 
on N, so we impose necessary restrictions on g G H, and 
determine a generalized coercivity condition that covers the 
case of the degeneracy. 



4.2. Existence, uniqueness, and wellposedness 

Before we investigate the existence and uniqueness of the 
solution a, we prove the following technical, but important, 
lemma. 

Lemma 4.1 Let P be an orthoprojector in H with a null- 
space N(P) = N and a range R(P) = P = N^, and A be 
a linear and bounded operator such that < A* = A on H. 
Then 

N(P ± A) n N = N(A) n N, (10) 

{N(P- L yl) n N} 1 - = R(A) + P, (11) 

N = (N(P- L A) n N) © P-LRL4). (12) 

Proof. We first verify (10). It follows from N(P ± A) D N{A) 
that the right-hand side of (10) is included in the left-hand 
side. To prove the reverse inclusion, let (p G N and P^Af = 
0, then = (P^Aip^ip) = (Atp,ip) = p^Vll 2 (recall that 
A > 0). Then A 1 / 2 : ip = and Aip = 0. Therefore, equality 
(10) holds. 

Equality (11) follows from (10), by substituting N(A) = 
F and R(A) = F in the well-known simple identity 
F i nP i = (F + P) x and noting that (R(A) + P) xl - = 

R(A) + P = R(A) + P by properties of the closure. 

Finally, to obtain the second term in the orthogonal de- 
composition _(12)_of_N we see that by (11) {N(P^A) n 
N}^ n N = R(A) + P n N; at the same time 

TL(A) + P n N = P-LR(A) + P n N 

= (p^rJa) © p) nN = P 7T RjAj, 

which completes the proof of the lemma. □ 

We start with the solution uniqueness. 
Lemma 4.2 Suppose that for some fixed g G H and f G H 
there exists a solution a of (8). Then it is unique provided 
that N(A) n N = {0}; otherwise, all possible solutions yield 
the hyperplane a + {N(A) n N} and there exists the unique 
normal (with minimal norm mHJ solution of (8) that can 
be also defined as a common element of the above hyperplane 
and the closed subspace R(-A) + P, which is the set of all 
normal solutions for all possible f and g. 
Proof. All solutions of (8) with g = f = constitute 
the closed subspace N(P ± A) n N(may be 0-dimensional), 
which by (10) is the same as N(A) n N. Hence, all solu- 
tions of (8) with the given g and /, provided that there 
exists at least one solution a, constitute the hyperplane 
a + N(A) n N. It is known that each closed hyperplane 
in a Hilbert space has a unique element with the minimal 
norm, i.e. the element that is orthogonal to the directing 
closed subspace N(^4) nN of the hyperplane. The orthogo- 
nal complement to the directing closed subspace is already 
given by (11). □ 

In the rest of the subsection we use the following equation 
equivalent to (9): 

(P ± A + P)a = P ± g + Pf. (13) 



r. 



The assumptions on the right-hand side of the system (8) 
which ensure the existence of a solution arc rather standard 
and follow from (13) easily. 

Lemma 4.3 For any f G H there exists a solution of (8) if 
and only if g G R(A)+P, i.e. P^g + Pf G P J -H(A) + P = 
R(A) + P. 

Proof. The subspace (not necessarily closed) P ± R(A) + P 
is simply the range of the operator P^A + P of equation 
(13). □ 

The subspace R(A) + P that appears in Lemmas 4.2 and 
4.3 plays the central role in the following necessary and 
sufficient conditions of wellposedncss. 
Theorem 4.1 The following statements are equivalent: 
(i) The subspace R(-A) + P is closed, 
(ii) The subspace AN + P is closed. 
(Hi) The subspace P- L H(A) is closed, 
(iv) The subspace P^AN is closed. 

(v) Problem (13) with f G H and g G R(A) + P is well- 
posed in the factor- space, ||cr||H/{N( J 4)nN} < c (ll.9ll + 
ll/H), or, equivalently, \\a\\ < c(\\g\\ + ||/||) for the 
normal solution a G R(A) + P- 
Proof. (1)^(3) We have R (A)+ P= P ± R(A) 8 P. 

(1) ^(2) The subspace P X K(A) P = K(A) + P is the 
range of the operator P^A + P. The range of a bounded 
operator is closed if and only if the range of the conjugate 
operator is closed. 

(2) o-(4) Using the same arguments as above, AN + P = 
P^AN 8 P. 

(1)<=>(5) The operator P^A + P is bounded and defined 
everywhere on a Hilbert space, thus it is closed. Therefore, 
the (pseudo)inverse operator 

(P ± A + Py 1 : K(P J -A + P)^ H/{N(A) n N} 

is closed. It is bounded if and only if its domain of definition 
R(P^A + P) is closed. A normal solution is a convenient 
representant of a factor-class in a Hilbert space. □ 

4.3. Generalized coercivity conditions 

Statements (l)-(4) in Theorem 4.1 may not be so easily 
verifiable in practice, so we want to find a somewhat eas- 
ier assumption that generalizes the standard coercivity as- 
sumption A > c a I on N C H, which itself does not hold if 
the operator A vanishes on a nontrivial subspace of N C H. 

Let us return back to equation (9). We remind the reader 
that the first equation in (8) is equivalent to the orthog- 
onal expansion a = ip + Pf, where ip = P^a G N. This 
and the second equation in (8) lead to (9) that we present 
here, introducing a special notation K = P^A |n, in the 
equivalent form 



Kip = <j>, ip G P- L H(A), <p = P ± g - P ± APf G P X TL(A).(U) 

under the assumption that g G R(A) + P. 

The operator K is bounded, selfadjoint, and nonnegativc 
definite on N, where N C H inherits the scalar product 
and the norm of H, so there exists a bounded, selfadjoint, 



and nonnegative definite square root \f~K on N. Applying 
the inf-sup condition to the operator \J~K on N, by direct 
analogy with Lemmas 3.1 and 3.2 and their proofs, we have 

that N {^/K) = N(if) and 

Theorem 4.2 The following statements are equivalent: 
(i) The subspace R (^/K^j C N is closed. 

(ii) The subspace R(A') C N is closed. 
(Hi) The inf-sup condition for the operator \f~K on N 
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(15) 













- = inf — 




N 



inf sup 

£SN CTeN ||e||N/N(if)||0-||N f £N 

>0 



N/N^) 



holds. 

(iv) The norm of the operator K^ 1 : R(A') — > N/N(A') 
is equal to p < oo. 
Moreover, under either of the assumptions we have 
r(Va) =R(A). 

Noticing that R(A) = P ± AN, we immediately see that 
statements (4) in Theorem 4.1 and (2) in Theorem 4.2 are 
the same, so all statements of Theorems 4.1 and 4.2 are 
equivalent. Our last goals in this subsection are to present 
statement (3) of Theorem 4.2 in original terms, so that it 
resembles the coercivity condition, and to bound the norm 
of the solution in terms of the norms of the right-hand sides, 
using statement (4) of Theorem 4.2. 

Theorem 4.3 For any g G R(A) + P the following as- 
sumption 



A > —I on the subspace P J "R(A) 
P 



(16) 



with a (finite) constant p > is necessary and sufficient for 
the normal solution a with P^a G P^R(A) to exist and 
to be unique and continuous in f G H and g G R(A) + P. 
Moreover, assumption (16) implies 



H 2 < ||/|| 2 + p 2 ||. 9 -AP/|| 



(17) 



Proof. First, we note that inequality (16) on the subspace 
P R(A) is equivalent to the same inequality on its clo- 
sure P J -H(A) because of the continuity of A and the scalar 
product. Second, as (e,Ke) = (e,P Ae) = (e, Ae) for all 
e G P ± R(A) C N, inequality (16) is also equivalent to 



K > —I on the closed subspace P ± R(A). 
P 



(18) 



Now we show that (18) is equivalent to (15), which is 
condition (3) of Theorem 4.2. For the numerator in (15), 

we have 



Ke = (e, Ae). To handle the denominator in 

N 

(15), we remind the reader the orthogonal decomposition 
N = (N(P ± A) n N) 8 P-LR(A) stated as (12) and proved 
in Lemma 4.1. Splitting e G N according to this orthogo- 
nal decomposition, we see that its first component — from 
N(P ± J 4) n N = N(A') — vanishes both in the numerator, 
since it is in the null-space of A, and in the denominator 



of (15), by the definition of the factor-norm, which gives 
(18), where only the second component — from P-'-R(A) — 
survives. 

We conclude that (16) is equivalent to (15), which is 
condition (3) of Theorem 4.2, and thus, to all statements 
of Theorems 4.1 and 4.2. Finally, if (16) holds then the 
subspace H(P ± A) is closed, the operator K : P ± R(A) i— > 
P- L H(A) is an isomorphism and problem (14) is wcllposcd 
for / G H and g € TL(A) + P, i.e. 



HVII < P\\P^9 - P ± APf\\ < p\\g - APf\\ 



(19) 



by Theorem 4.2. Estimate (17) follows from a = ip + Pf 
and (19) due to the statement of Lemma 4.2 that the nor- 
ma l solution a E R(A) + P, that is, ip e P ± H(A) = 
P ± Ti(A) + P = P ± n (P 1 - n N(A))- 1 is the corresponding 
part of the orthogonal expansion a = ip + Pf for the nor- 
mal solution. □ 

4.4. Minimum gap between subspaces 

The rest of the section concerns the case where the range 
of A is closed, so assumption (16) can be equivalently re- 
formulated using the minimum gap between some relevant 
subspaces. We first find a simple way to check if the range 
of A is closed. 
Lemma 4.4 Condition 



A > — I on the subspaceR(A) = D 

PD 



(20) 



with a (finite) constant po > is equivalent to closedness 
o/D. 

Proof. The operator A is a linear, bounded, and every- 
where defined. Thus, it is closed and its inverse A -1 : D — > 
H/N(A) is also closed. Boundcdncss of the inverse is equiv- 
alent, on the one hand, to condition (20) and, on the other 
hand, to closedness of D. □ 

Now we are ready to present a simplified version of the 
necessary and sufficient condition of wellposcdness (16), 
assuming that the range of A is closed. 
Theorem 4.4 Let the range R(-A) = D be closed, the or- 
thoprojector on D be denoted by D, and the constant pu > 
be defined by (20). Then inequality (16) is equivalent to the 
inequality 

K = inf W>0- (21) 

In particular, (20) and (21) lead to (16) with p = pd/k 2 - 
Proof. We have R(P ± A) = F- L R(A) = P X D, i.e. the 
subspaces indicated in (16) and (21) coincide. Now the 
main assertion of the Lemma is a consequence of relations 
J-(DtJj,D^) < (Atp,i/)) = (AD?P,DiP) < \\A\\(Dil),D$) 
which hold for an arbitrary tp G H. □ 

The next two lemmas provide alternative assumptions, 
equivalent to (21), which are necessary and sufficient for 
wellposcdness, assuming that the range of A is closed. It is 
important to have a choice of a criterion that may be easier 



to check in a practical application. For aesthetic reasons we 
denote N = P x . 

Lemma 4.5 Let D and P be orthogonal projectors onto 
closed subspaces D and P, and let D = I — D and P x = 
I — P be orthogonal projectors onto the orthogonal comple- 
ments D x and P , respectively. The following statements 
are equivalent: 
(i) The subspace P X D is closed. 

(ii) The subspace D + P is closed. 

(Hi) The subspace D 1 - + P x is closed. 

(iv) The subspace PD X is closed. 
Proof. (1)<=>(2) The subspace P ± D is closed iff the sub- 
space P ±W D © P = D + P is closed as the terms are orthog- 
onal in the first expression. 

(2) ^(3) By Theorem IV-4.8 of [12], a sum of closed sub- 
spaces in a Hilbcrt space is closed if and only if the sum of 
their orthogonal complements is closed. 

(3) o-(4) Using the same arguments as above, P x + D x = 
P x 8 PD X . □ 

Lemma 4.6 Using the notation of Lemma ^.5, the follow- 
ing equalities hold: 



inf 



dist{?/>; D} 



inf 



<tist{ip; P} 



V>eP, V£D dist{?/>; DnP} V£D, </>£p dist{V>; PnD} 

. nf dist{^; P^} 

i/>eD^, ^p^ dist{-0; P ± n D^} 
. f dist-Q; D x } 

_ y>eP-<-, V'^d^ dist{?/>; D x n P x } 



= inf 



\Lhl)\ 



inf 



D^ip\\ 



inf 

'ipeD ± p 



\pn 



= inf 



Moreover, each statement in the previous Lemma is equiv- 
alent to the positiveness k > in (21). 
Proof. The first three equalities are derived in Section IV-4 
of [12] on the minimum gap between subspaces, along with 
a statement that positiveness of the minimum gap between 
two given subspaces is a necessary and sufficient condition 
of the sum of the subspaces, in our case, D+P, to be closed. 
We now prove that 



inf 



dist-^D^} 



■(igp 1 , ^d 1 dist{V>; D- 1 - n P- 1 } 



ipeP ± r>\{o} \\ip\\ 



All other equalities can be then trivially derived from the 
previous ones just by interchanging P and D. 

Wc first notice that in the right-hand side we can apply 
the inf to the closure P^D \ {0} as well, because a norm is 
a continuous function, 



inf 

V-eP J -D\{o} 



\Dj'\ 



inf 

t/ieP 1 D\{o) 



\D4'\ 



We have, P X D = P x n (P 1 - n D 1 ) 1 as N(Z?P- L ) = P 
(P x n D- 1 ). The latter can be checked directly. 



<7 



Wc always have dist{V>;D^} = \\Di(>\\. If tp G P^D = 
P- L n(P- L nD ± )- L C (P- L nD ± ) ± ,wc also have distjV'iD-Ln 
P- L } = ||^||.Thus, 

dist-j^D- 1 } _ \\Dij}\ 
distil; D-LnP 1 } ~~ ij^jj 



^ G P^D^O}. 



Finally, using the orthogonal representation P 1 - = (P n 
D J -)©P- L D, every <p G P- 1 can be written as the orthogonal 
sum 99 = (95 — ip) ® V: where (p — ip£ P- 1 P\ T) 1 -, ip E P-"-D. 
Then dist{^; D^} = dist{<y3; D- 1 } and also distj^D- 1 n 
P- 1 } = dist{(^; n P- 1 }; so the value of the ratio 



dist^D^} 
dist{V>; D 1 n P 



dist{ip;D^} 
^ ~ dist^; D-LRP-L} 



does not depend on ip — ip and its two infimum values, taken 
with respect to ip G P^D \ {0} and ip G P ± , <p <£ D x , 
coincide. □ 

Finally, we notice that g = if we apply a saddle point 
approach to diffusion or linear elasticity equations. Indeed, 
in the Hellinger-Reissner formulation of nonhomogeneous 
Lame equations, our a represents the stress tensor, the La- 
grange multiplier u is the displacement, and if we also in- 
troduce the stain e by the stain-displacement relation e = 
—B*u, then the first line in system (5) becomes Aa — e = 
g, which is the constitutive equation (3-D Hooke's law), 
where of course g = 0. The second line in (5) is the equi- 
librium equation, where all body and traction forces are 
represented by / 7^ 0. The assumption 5 = allows us to 
look for even weaker conditions of wellposedness that we 
plan to investigate in the future. 
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